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ABSTRACT 


A modal solution of the flutter and vibration problems for a multiple 
component system is presented. The formulation requires the aerodynamic 
characteristics in the form of aerodynamic influence coefficients, and the 
structural and inertial characteristics in the form of free vibration modes 
and frequencies and a mass matrix for each component. The use of a rigid- 
body modal matrix permits a general analysis for a system free in space 
with up to six rigid-body degrees of freedom. The solution also utilizes a 
newly developed subroutine for finding the eigenvalues and eigenvectors of a 
non-Hermitian matrix having close roots, and obtains the corresponding 
values of frequency, required structural damping, and velocity of the flutter 


system. 


The Aerospace IBM 7090 Computer Program No. LDO14A provides 
the flutter or vibration solution for a system composed of as many as 20 
flexible components, eachofwhich may have upto 50 controlpoints. The solu- 
tion is carried out by a modal (series) method choosing a series of up to 
40 cantilever and/or free-free vibration modes to approximate the control 
point deflections. The program may be used to find the vibration character- 
istics of a composite system if the cantilever vibration characteristics of 
the components are known, and the generalized masses may be obtained for 
use in subsequent modal analyses of flying qualities or stability and control 
of flexible vehicles. An option is provided to vary the density as well as the 


reduced velocity in the flutter analysis. 
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SYMBOLS 
A Element of aerodynamic force matrix, lb 
a Element of normalized flexibility matrix, in/1lb 
a Amplitude coefficient of cantilever free vibration modal series, in. 
a Amplitude coefficient of generalized centroidal coordinates 
ap Generalized amplitude coefficient of rigid-body modal series 
ang Element of rigid component modal matrix 
b. Reference semichord, ft 
Ch Element of oscillatory aerodynamic influence coefficient matrix 
F Number of control points at which the modes are specified 
g Structural damping coefficient 
ho Control point deflection due to rigid-body motion, in. 
hy Control point deflection, in. 
K Flexibility matrix normalizing constant 
Kp Element of generalized spring matrix 
k Element of cantilever spring matrix, lb/in. 
k Reference reduced frequency 
M Element of mass matrix, 1b 
M Element of complex mass matrix including aerodynamic effects, lb 
m Element of generalized mass matrix 
™m Element of sum of generalized mass and aerodynamic matrices 


Number of assumed flexible modes 
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SYMBOLS (Continued) 


Element of generalized aerodynamic force matrix 
Number of rigid- body degrees of freedom 
Reference semispan, ft (i.e., span measured from root to tip) 
Velocity, knots 

Element of aerodynamic weighting matrix 
Cartesian coordinates 

Error of series solution 

Air density, slopes? 

Eigenvalue, (opel ° 

Frequency, cps 

Prime denotes use of free-free modes 

Square matrix 

Column matrix 

Diagonal matrix 

Transposed matrix 


Inverse matrix 
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SECTION I 


FORMULATION OF PROBLEM 


A, Introduction 


The use of structural influence coefficients in a collocation method of 
vibration analysis and of both structural and aerodynamic influence coe ffi- 
cients in a collocation flutter analysis (both of which are reviewed in Ref. 1) 
provides a number of advantages in terms of accuracy and convenience. 
However, the computation time in handling very large matrices becomes 
prohibitive, while computer capacity also may limit the analysis of extremely 
large systems, and a series or modal solution must be used. A modal for- 
mulation of the flutter problem using aerodynamic influence coefficients has 
been presented in Ref. 2. The present report reviews the derivations in 
Ref. 2 and extends its computational aspects, particularly in regard to the 
simplifications that accompany an option for a modal vibration analysis and 
an option to vary the density as well as the reduced velocity in the flutter 
analysis. The present program also utilizes a recently developed subroutine 
for finding the eigenvalues and eigenvectors of a non-Hermitian matrix having 
close roots, and its requirement for modes and frequencies is consistent with 


those generated by the vibration option of Ref, 1. 


B. Sign Convention 

The NASA body axis system with the x, y, and z axes directed forward, 
starboard, and downward, respectively, is recommended for consistency with 
the formulation of the static aeroelastic problems in Ref. 3. However, the 
usual flutter convention with the x, y, and z axes directed aft, starboard, and 
downward, respectively, may be used instead. In either case, the vertical 


normal force and deflection are positive downward. 
Cc. Derivation of Equations 


It is convenient to present two derivations of the modal formulation of 


the flutter and vibraticn problems according to whether cantilever or 


free-free modes are to be utilized. The design analysis normally will make 
use of cantilever modes, but often the final analysis will be able to take ad- 
vantage of the free-free modes available from ground vibration measurements. 
The two derivations are only slightly different; nevertheless it is useful to 


separate them. 


The modal formulation may be derived equally well from an assumed 
series for the deflection modes and either Lagrange's equations for the 
dynamic system or the collocation (lumped parameter) equations of motion 
and Galerkin's method. Since the collocation equations of motion have been 
developed in Ref. 1, it is convenient to choose the latter approach. It is 
further convenient to summarize the significant features in the derivation of 
Ref. 1 to provide the basis for the present derivation. In Ref. 1, a single 
type of coordinate, viz. , the deflection h, has been chosen as an adequate 
measure of both the deformation and the free-stream disturbance, not only 
for simplicity in the resulting equations but also because deflections have a 
more general meaning on a cambered vehicle and deflection influence coeffi- 
cients are more readily obtained from a structural analysis than slope (or 
twist) influence coefficients. The resulting deformation matrix equation 


becomes 
{ay} - fag} = KALE} + FA) q) 


where {hy | is the set of components of the absolute deflections of the control 
points, {ho} is the set of components of the deflections of the control points 
due to the rigid-body motion of some reference point, [a] is the set of struc- 
tural influence coefficients (SICs, or flexibility matrix) for the system 
cantilevered from (or otherwise restrained at) the reference point, {F,| is the 
set of inertial force components integrated throughout the region adjacent to 
each control point, IF} is the set of aerodynamic force components integrated 
over the vehicle surface adjacent to each control point, and the scaler K has 
been introduced as a factor to the SICs for convenience in investigating varia- 


tions in stiffness levels. 


The inertial forces may be written in terms of a mass matrix [M] and 


the control point accelerations. 
{Fr} = -(1/386)[M]{h, } (2) 


where the diagonal elements of the mass matrix are found from integrating 
the structural mass density throughout the region adjacent to the control 
points. (N.B. The mass matrix need not be diagonal, and, in general, will 
not be so if the elements must be derived from a set of weight data previously 
lumped at a system of control points different from those required in the 
aeroelastic analysis. The use of a coupled mass matrix permits simulation 
of given inertial characteristics at a set of control points frequently dictated 
by more difficult aerodynamic considerations. A special technique for 
deriving the elements of a coupled mass matrix corresponding toa specified 
set of control points from a discrete set of given weight data is suggested in 
Appendix B of Ref. 1.) 


The aerodynamic forces may be written in terms of a set of aerodynamic 
influence coefficients (AICs) for oscillatory motion and the control point 


deflections as 
[F,] = (407/12) pubs (WILC,}{hy} a 


where [c,] is the theoretically derived dimensionless (complex) matrix of 
oscillatory AICs, wis the frequency of the assumed harmonic motion, p is the 
atmospheric density, b, is the reference semichord, s is the reference span, 
and [W] is an empirically derived weighting matrix for modification of the 
theoretical AICs. A method for obtaining the elements of the weighting 


matrix has been suggested in Ref. 4. 


The sum of the force components may now be written from Eqs. (2) and 


(3) for the case of harmonic motion. 


{Fi} +{F,} = (407 /386)([M] + 32.174 pb“s[ W] [C.]){h,} (4a) 
= (4nw°/386)([M] +[A]){h,| (4b) 
= (4nw"/386)[M] {h,| (4c) 


We next discuss the manner of inclusion of the rigid-body degrees of 
freedom in Eq. (1). The matrix {ho} has been defined as the set of compo- 
nents of the deflections of the control points due to the rigid-body motion of 
the reference point. Each component of the control point deflections ho is 
linearly related to the rigid-body translations and rotations provided the 
rotations are small. Therefore, we may define a rigid-body modal matrix 


[hp] as the transformation 


{ho} = (hp lag] (5) 


where {ap} is the set of amplitudes of rigid-body translations and rotations of 
the reference point. As an example, if we consider symmetrical vertical 
motion, [hp] is composed of two columns: The first is a unit column corre- 
sponding to the plunging mode; the second consists of the x-coordinate of each 
control point corresponding to the pitching mode; {ap | is composed of two 
elements: The first is the plunging displacement Zz) the second is the pitching 


angular displacement @. 


Before proceeding in the derivation, we should review the format of the 
various matrices in the case ofa multiple flexible component system. As an 
example, we consider a symmetrical flutter analysis of an aircraft whose 
wing, aft fuselage, and tail are flexible and whose forward fuselage may be 


assumed to be rigid. We assume that the reference point (the cantilever 


a a 


point) can be located in the vehicle such that its various components are 
independent. If we choose a point at the intersection of the wing and fuselage, 
then the wing is independent of the aft fuselage-tail combination, but the tail 
and aft fuselage must be considered together. The motion of the rigid 
forward fuselage is determined by the motion of the reference point, and the 
forward fuselage will not enter into any flexible considerations but only into 
the free-free boundary conditions. From the foregoing, it is seen that the 
various matrices will appear in partitioned form. If we denote the wing and 
aft fuselage-tail system by the subscripts 1 and 2, respectively, then the 


flexibility matrix appears as 


ay 0 
fa] Eee osPeeoe (6) 
0 as 
the mass matrix as 
M, 0 
[M] = |----- ease A) 
1 
0 M, 
the weighting matrix as 
Wy 0 
oo ia chad aaa (8) 
1 
¢) WwW, 
the AICs as 
Chi 0) 
[c,] = ‘sepeaea ----- (9) 
0 Ch2 
a5 


and the rigid body modal matrix as 


Pea 


[hp] cag | (aes (10) 


hpe 


Two requirements should be emphasized with regard to the AlCs. The first 


concerns the proper inclusion of the reference geometry associated with the 


nondimensional AICs. The dimensional form of Eq. (9) may be written 


z 
brs{C,] Sle SSS se peg eters? (11) 


where b,. and s are the reference semichord and span of the composite 
system, b, and $s, are the reference geometry for the first component, and 
by and S, are the reference geometry for the second component. The second 
requirement is that the AICs for each component must be determined for the 
same "dimensional" reduced velocity V/w. If the reference reduced velocity 
is 


1/k_ = V/b_w (12) 
r r 


then the reduced velocity for the first component must be 


1/k, = (1/k Xb, /b)) (13) 


and for the second component 


1/k, = (1/k, ib _/b5) (14) 


=e 


Both of these requirements can be met in formulating the composite AICs by 
choosing the same reference geometry in determining the AICs for each 


component. 


The rigid-body modal matrix provides the basis for a general statement 
of the boundary conditions for the free-free flutter of the composite system. 


The boundary conditions for harmonic motion may be written as 
[hy]*[M]{h,} + [a7] {ap} =o] (15) 
R 1 R 


where [Am] is an incremental generalized mass matrix, including aerodynamic 
effects, of any rigid component of the system attached to the reference point 
(e.g., the forward fuselage that was assumed to be rigid in the foregoing 
example), and not considered in the formulation of the flexible component 
mass and aerodynamic matrices. The form of the matrix [Am] may be illus- 
trated by the previous example with the rigid forward fuselage again in sym- 


metrical motion. We may write 
[am] = [Am] + [AQ] (16) 
where the generalized rigid component mass matrix of the forward fuselage is 


M NS) 
fo} (eo) 
[Am] = (17) 
Ss 


I 
fe) yo 


in which M,: So: and I 56 are the mass, static unbalance about the reference 
point, and pitching moment of inertia about the reference point, respectively, 
of the forward fuselage, and the generalized aerodynamic forces on the 


forward fuselage (if not negligible) are found from 


[AQ] = 32.174 pb ath, (eG. [np] (18) 


where [hp] is the rigid-body modal matrix, [CL] is the set of AICs, and Be 
and S are the reference geometry for the forward fuselage. Again the AICs 


must be found for the reduced velocity of the composite system. 


The collocation equation for flutter follows by substituting Eq. (4c) into 
Eq. (1) and adding the structural damping factor 1/(1 + ig) to the flexibility 
matrix to provide the artificial structural damping necessary to sustain the 


assumed harmonic motion of the flutter system 
2,- 2,: : <G 
{h} - {ho} = (40 Ko"/386 (1 + ig)al{M}{h,} (19) 


The rigid- body degrees of freedom may be eliminated by use of Eqs. (5) and 
(15). For the purpose of seeking series solutions to Eq. (19), it is conven- 


ient to introduce a cantilever spring matrix defined by 


(k] = (386/K)[a]"} (20) 


and to rewrite the flutter equation as 
27 
o{k]({h,} - [h}) = 47 [M]{h,} (21) 


where the eigenvalue 2 = (1 + er. 


This review of the derivation of the collocation equation and of the 
formats of the various matrices places us ina position to consider the two 
modal formulations of the flutter and vibration problems. We consider first 
the choice of a series of cantilever modes, and then the choice of a series of 


free-free modes. 


Jee. 


Case I: Cantilever Modes 


Let us denote by N the number of flexible modes given, by F the number 
of control points at which the modes are specified, and by R the number of 
rigid-body degrees of freedom. If we select a series of cantilever modes 
(which satisfy the requirement of the Galerkin method that all boundary con- 


ditions be met) for the solution of Eq. (21) then the series may be written 
{hy} - {ho} = (hp lay} (22) 


where [nh] is the matrix of given cantilever modes and the {arn} are the modal 
amplitudes. The matrix [h J appears in a partitioned form (size F X N) in 


which each partition corresponds to a component of the composite system 


[hp] = eae (23) 

- ! hy? 

Substituting Eqs. (5) and (22) into Eq. (15) yields 
[mp elliapl + [Mprllap!} = 0 (24) 

where 

[ppl = [mppl + (Qggl (25a) 
= (hp) [Mth] + fhp} [Althg (25b) 
[Tiga] = [mpg] + [Opp (26a) 
= (fnpi [MIlag] + f4m]) + ([hpFTAlthgl + (40) (268) 


and from which 
lap} = -[figgl [gpl p! (27) 
The flutter equation, Eq. (21), now becomes 
Ok][h play} = 40MM (a pl{ap} + Cagllag}) + te] (28) 


where {e} is the error in the series solution. Applying the Galerkin principle 
(Ih l fe} = 0) yields : 


AK, lap) = 40 ppllap} + pRl{ap)) (29) 
where 
[K,] = (hy) [kihy] (30) 
[mop =[™ppl + [Qe (31a) 
= (hb, tM) (ny) +[hp) “Ll fay) (31) 
(Mop) =[mpp] +[QpR] (32a) 
=[mp pl? +{hpl LA] {hg (32b) 


216:- 


” 


If we eliminate the rigid-body degrees of freedom by substituting Eq. (27) 


into Eq. (29), we obtain the final form of the cantilever modal formulation 
UK aap} = 40 (mMppl - Hppllpel (mp -)lap| 33 
Spier! > FF FRIL™RR RF) Ip (329 
or in the canonical form 
ofan} = (K,/407) iippl - FppliFiggl gpllagl (34) 
a a FF FRO RR RFR 


We note that the generalized spring matrix is diagonal because our 
formulation is not statically coupled (it is dynamically coupled instead). The 
generalized spring matrix may be found in terms of the generalized free 


vibration mass matrix and the natural (normal or uncoupled) frequencies 
[K [an] = [diagonal elements of [m..]] ioe) (35) 
F 8 FFI °F 


The flutter equation, Eq. (34), can be solved by the power method for 


the eigenvalues, 2 = QR + iQ, from which we find the frequency 


w = NOR (cps) (36) 
the required structural damping 
= 2,/Qp (37) 


and, since the aerodynamic matrix required the assumption of a reduced 


velocity, 1/k,, for its formulation, the velocity 


V=0. 5921(2Twb_) 1/k,) (knots) (38) 


a 


The power method as developed for non-Hermitian matrices with close roots 
is discussed in Appendix A of Ref. 1. The flutter mode is found from the 


originally assumed series 


{hy} 5 [(rLliap} a [np ]{apg} (39a) 
= ({hpl - (bp Figg) Tpelayl (39D) 


From a series of solutions of Eq. (34) for various Pedweed velocities 
the conventional required damping versus velocity stability curve can be con- 
structed for a specific altitude and the flutter point is determined as the 
velocity for which the required damping and actual damping are equal. An 
alternative approach to the flutter analysis is based on a single representative 
reduced velocity and a series of solutions of Eq. (34), carried out for various 
densities. The density at which the required damping and actual damping are 
equal may be used to find a stiffne ss-altitude similarity parameter for flutter 
from which the flutter stability may be determined. However, at present, the 
validity of this latter approach requires further investigation, particularly 


the sensitivity of the results to the choice of repre sentative reduced velocity. 


The generalized mass corresponding to each free vibration mode is also 
of interest in various modal analyses of flying qualities, stability and control 
characteristics, or transient response of flexible vehicles. If Eq. (34) is 
solved for the free vibration modes (by deleting the aerodynamic terms), then 


the n'th generalized mass is given by 


it). {nl /7 aya) + fal?) 7 fam]fal?)) (40a) 


N.B. When the vibration modes {ni} are normalized on their maximum 
elements, the modal amplitudes jal} cannot be normalized but must be 


scaled accordingly for the calculation of the generalized masses. 


ee 


mi") = fal) Fimppllae| + lag i imppllan | 
n aS i = (40b) 
+ fal) Pimp pale) + (ae) T img ella 
= {a ay i [mpF] (ale?) i [mp pl a) (40c) 


where fa y is the n'th set of cantilever modal amplitudes and the corres- 
ponding rigid component mode is found from Eq. (27) (again without the 


aerodynamic terms). 


=-[mpal ‘fp play | (41) 


Case II: Free-Free Modes 


When free-free vibration data are available for a flutter analysis, itis 
necessary to write two series expressions that provide the basis for the 


modal solution. 
{hy} = (hh lah} + fhgllag] (42) 
{ap} = lappllapl + {ag} (43) 


where the prime denotes the free-free condition, {aq} is the matrix of 
centroidal generalized coordinates, and [ape is the rigid component modal 
matrix. The rigid component modal matrix may be found from the boundary 
conditions for the free vibration. However, it is desirable to permit the 
consideration of some cantilever modes in addition to the orthogonal free-free 


modes. Hence, if we write 


Be 


eae. Pht] (44) 
then 
famed = foe. 10 (45) 


where the subscripts o and c denote orthogonal and cantilever, respectively. 
The orthogonal partition follows from the free vibration boundary condition 


and is given by 
] = -[4m]" "fag EMI] (46) 


[an‘o 


Making the series substitutions, Eqs. (42) and (43), into the boundary condi- 
tion, Eq. (15), yields 


[hp IMM Chat} + fhgllag}) + map gllat} + {ag} = (47) 
from which 
fag} = gal [Rell ete) 
where 
(mt, -] = (hg) [Mb] + 4M ]lap yl (49) 


In order to apply the Galerkin principle when free-free modes are used in the 


series solution, it is necessary to obtain the equation of motion of the rigid 


44 


tM my 


||! Cae F 2 


component. If the flexible system mass term is eliminated between Eqs. (15) 


and (21), then the rigid component equation of motion is obtained. 
-Ah,} [kl{h,} - {h_}) = 407[amlia, | (50) 
R uw o R 
Substituting the series solutions, Eqs. (42) and (43), into Eq. (50) yields 
Ah, F (kh at.} = 40am (ap lata} + {ac} + fen} (51) 
R Facer RFVIF G R 
where ler] is the error in Eq. (50) due to the series solution and 


[ne] = (hi) - (52) 


[hp laps] 


Substituting the series solutions into the flexible system equation of motion, 
Eq. (21), yields 


Akh" at} = 4n7(M] Chat} + (haWag}) + fes| (53) 


where {em his the error in Eq. (21) due to the series solution. Applying the 
Galerkin principle (in this case, laeel ler] + [h! ml Te } = = 0) yields the 


generalized equation of motion 
a = 
AKER] fap] = 47° U®ppllap| + pgllag) (54) 
where 


[Kip] = (a) Cxl[hig] (55) 


Ape 


Mp] = (mppl + (Qp el) (56a) 


= ({ht J (MI[hL] + [appl mag pl) 
+ (ht) CADE] + agp) Oleg pl) (56b) 
[Mito] = (ht (Mlthp] + [appl [am (57) 


If we eliminate the centroidal degrees of freedom by substituting Eq. (48) 


‘into Eq. (54), we obtain the final form of the free-free modal formation 
; 2732) —_ eek -l-_ 
Oo Kata} = 40 (mie) - mpglMpggl [Rpllay| (58) 
or in the canonical form 
ofat.} =f Ki /40? | (Gmtip) - Fpl iigg) Rely (59) 
Fi F FF FR RR RF F 


We again note that the diagonal generalized spring matrix is found from the 


generalized mass matrix and natural frequencies 
[K'/4n2] = [diagonal elements of [m',p] Ifo!’ (60) 
F 8 FFINOp 


The solution of the flutter equation by iteration has been discussed in the 
previous case. The flutter mode follows from the assumed series. 


Substituting Eq. (48) into Eq. (42) yields 


(nh } = Ub) - (hg liigg) TH aellatp (61) 


sie 


he ) 


The n'th generalized mass follows from the solution of Eq. (59) with the 
aerodynamic terms deleted, and Eqs. (40a), (42), and (43). 


ml”) = fal) Pimt ela”) + lag? lag llag” 


+ {al Tim, Ham) + (a) Fim pllak™ ‘hs 
= ee) | + [mip] {ae} (62a) 
where 
(aD) = Impl “lm ella] (63) 
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SECTION II 


GENERAL DESCRIPTION OF INPUT 


A. Units 


The dimensional data required for each component consist of the mass 
matrix in pounds, the free vibration frequencies in cycles per second, and 
the reference semichords and semispans in feet. The free vibration modes 
and the aerodynamic influence coefficients are dimensionless. Compatibility 
is required between the units of the rigid-body modal matrices and the 
generalized rigid-body mass matrix; for example, if So and I__sare in the 
pound-foot system, then the x-coordinates, which correspond to the pitching 
mode, must be measured in feet, whereas if the generalized masses are 
given in the pound-inch system, the x-coordinates must be measured in 


inches. The density is required in slugs per cubic foot. 
B. Classes of Data and Example Problems 


Five classes of data must be provided: mass, aerodynamic influence 
coefficients (AICs) and their weighting matrices, free vibration modes 
(cantilever and/or free-free) and corresponding frequencies, rigid-body 
modes, and the generalized masses for any rigid component in the system. 
(The cantilever case does not require the rigid-body modes or the generalized 


masses). 


Two example problems will be considered. Each problem will be a 
symmetrical flutter analysis of the jet transport wing treated as an example 
throughout Ref. 5. In the first analysis, we will use cantilever vibration 
modes and frequencies, and in the second we use free-free symmetric modes 
and frequencies such as might be obtained from a ground vibration survey. 

In each problem, we assume that four modes will be sufficient to obtain 


convergence in the flutter solution. 
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The four flexible cantilever modes (calculated by the collocation method 


of Ref. 1) are 


[he] = 

0.045185 -0. 163556 -0. 219686 0.291415 
0.042502 0.544990 -0. 243887 0.279638 
0.148256 -0. 279631 -0. 205463 0. 306804 
0.143727 1.0 -0. 362884 0.353608 
0.379411 -0.226740 -0. 371824 -0. 559460 
0. 376733 0. 946563 -0. 428613 -0.490161 
0.692995 -0. 138889 0. 203824 -0. 146405 
0.692655 0. 887880 0. 163356 -0.020285 
0.998969 -0.039737 1.0 : 0.837729 
1.0 0.822370 0.981915 1.0 


The corresponding frequencies are 


, = 2.036802 cps 
w> = 3.552589 cps 
w3 = 7.280460 cps 
w4 = 11. 698543 cps 


The four free-free symmetric modes (also found by the method of Ref. 1) are 


[hp] = 
-0.129721 -0.091809 -0.035610 0.305979 
-0. 131332 0.568102 -0.047615 0.287803 
-0.017537 -0.221627 -0. 162942 0.442843 
-0.015998 1.0 -0.189775 0.503783 
0. 247806 -0.215371 -0. 449387 -0. 561462 
0.250893 0.905541 -0. 477740 -0. 443473 
0.622665 -0.205131 0. 106948 -0. 225098 
0.627940 0.776147 0.097215 -0.006219 
0.993742 -0. 186900 0. 984408 0.741022 
1.0 0.636687 1.0 1.0 
and the corresponding frequencies are 

w, = 2.442962 

@2 = 3.603115 

w3 = 8.553045 

W4 = 12.707073 
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The symmetrical rigid component (one-half of the fuselage) generalized 


mass matrix is given below in the pound-inch system. 


M Ss 17,400 1,370, 250 
fo) fe) 


[Am] = = 
S 


I 1, 370, 250 4, 457, 907, 200 
fo) yo 


(In the antisymmetric case, the rigid component generalized mass matrix 
is [Am] = [tls in the composite symmetric-antisymmetric case, the 


generalized mass matrix is 


M_._ Ss 0 
[e) ie) 
[4m] =] S, ore ) ) 
0) ) I =a 
xO 


The mass and aerodynamic matrices for the wing and the aerodynamic 
matrix for the fuselage are shown in the example problem printed output 
(pages 47-73) and will notbe repeated here. The aerodynamic weighting 


matrix is taken as unity, and no input is required. 


In addition to the flexible modes, the symmetrical analysis requires 
the two rigid-body modes for the degrees of freedom of plunging and pitching. 
The rigid-body modal matrix therefore consists of two columns: a unit 
column and a column of the x-coordinate of each control point. The rigid- 


body modal matrix for the wing is 


[ny] =[1 x] = 


a ee ol 
1 
fo 
Ww 
N 
[o) 


22) - 


The rigid-body modal matrix for the fuselage is 


1 -373.30 
1 -248. 30 
[h,! ={1 x] 1 -123. 30 
1 +1. 70 


1 +126.70 


Note that the above matrix is used in computing the incremental complex 
generalized mass which results from the aerodynamic loads on one-half of 
the fuselage and, therefore, the x-coordinates must be given in the proper 
order for the control points used in computing the fuselage AICs. (In this 
problem the fuselage AICs are hypothetical but the x-coordinates are given 
in the order necessary for the slender-body theory AICs referred to later 


(page 39) in this report.) 


The reference geometry for the wing is Dig = 5.468 ft and Sui = 37.917 ft. 
The reference geometry for the fuselage is Rig = 5.468 ft and le 18.9585 ft. 
(We have assumed that the wing reference geometry was used to nondimen- 
sionalize the fuselage AIC matrix, and since we need only one -half of the 
fuselage aerodynamic forces they can be obtained by setting = (1/2)s..-) 
The flutter analyses are carried out for the reduced velocity 1/k = 16.67 
using the reference semichord for the system bo = 5.468 ft. Both analyses 


are carried out for sea level; i.e., p = 0.002378 slugs per cubic foot. 


The results for each of the two problems are shown in the example 
printed output, Section IV, Part A. The final results are shown also in the 
cases of 6-, 8-, and 10-mode solutions to demonstrate the convergence of the 
four-mode solution. We note that the four-mode convergence is very poor, 
emphasizing the importance of checking convergence in any modal solution. 
We also note that the 10-mode solution agrees exactly with the collocation 


solution shown in Ref. 1 as is to be expected for a system with 10 degrees 
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of freedom. A final point of interest is to note that the use of free-free input 


modes does not provide any better convergence than the cantilever modes. 
on Program Restrictions and Options 


1. The maximum number of vibration modes that can be input for 


the system is 40. 


2. The maximum number of cantilever vibration modes is 40 


whether the system consists of one or more components. 


3. The maximum number of free-free vibration modes for the 
system is 40 minus the total number of cantilever modes to be used for the 
system. Since free-free modes represent the entire flexible system, a 


portion of each free-free mode must be input corresponding to each flexible 


component. 
4. The number of flexible components is S$ 20. 
5. The number of values used in the reference reduced velocity 


(1/k,) series is S$ 20. 
6. The maximum number of values used in a density series is £ 20. 


7. The program provides for varying the densities with each 1 /kk 


or for using the same densities with all 1/k 's. 


8. The maximum number of output modes is 25; however, the 


number cannot exceed the number of input vibration modes. 


9: The maximum number of control points for the rigid component 


or any flexible surface is 49. 
10. The number of rigid-body modes is € 6. 


ll. It is possible to reserve a partition in the upper left-hand corner 
of any flexible component AIC matrix for control points whose aerodynamic 
forces may be neglected or found from some alternate theory to that used for 
the primary control points. This partition is termed the "external stores" 


region since external stores are an example of a source of additional control 
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points requiring such special consideration. The maximum number of control 


points that can be reserved on each surface for external stores is S 48. 


12. A weighting matrix is optional input for each surface. The order 
of this matrix must agree with the order of the AIC matrices for the 


particular surface. 


13. Any number of complete sets (decks) of input data may be stacked 


and run in one machine pass. 
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SECTION III 


DATA DECK SETUP 


A. Loading Order 


Input data decks are loaded behind column binary deck No. LDOI4A. 


The data decks are assembled using cards punched from keypunch forms 


and cards that are punched-card output from appropriate IBM programs. 


The data items in each deck have the following order with this exception: 


some of the items may be omitted if indicators used in the control cards 


(Items 2 and 4) specify their absence. 


1, 
2. 
3. 
4, 


5: 
(1/k's). 


6. 


Title card. 

Data deck general control card. 

Data card for changes in matrix iteration tests. 

Control card(s) for external stores and weighting matrices. 


Reference semichord (b,) and reference reduced velocities 


Reference semichords (b 3) and reference semispans (s;) for 


rigid and flexible components (surfaces). 


7. 
8. 
9. 
10. 
a 


12. 


Density series (if same densities are used for all 1/k,'s). 
Generalized mass matrix ([ Am]) for rigid component. 
Mass matrix ([ MJ.) for each flexible surface. 

Rigid-body modal matrix (hp o)) for rigid component. 
Rigid-body modal matrix ([hp];) for each flexible surface. 


Frequencies ,) and vibration modal matrix (hp) i) for each 


flexible surface. 
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Weighting matrices ([ w],) : 


a25 = 


14. Aerodynamic input repeated for each 1/k., 


a. . Density series cards (if densities vary with each 1/k.), 
b. AIC matrix ([ Cho!) for rigid component (if present), 
es AIC matrix ([c,],) for each flexible surface. 
B. Input Data Description 
ae The title card may contain any alphanumeric characters 


desired in Columns 2 through 80. Column | should be blank. 


oe Data deck general control card (FORMAT 1814): The first 11 


fields in this card contain, in the order defined, the following indicators: 
NSUR: Number of flexible components (surfaces), £20. 


NAERO: Number of reference reduced velocities, £20; 
NAERO = 0 is used for the vibration analysis. 


NRIGID: Number of rigid-body motion modes to be input (= number 
of columns in [hp ]); NRIGID = 0 for the cantilever case. 


NFUS: NFUS must = 1 if AICs (oo): j = 1, NAERO) are input 
for the rigid component. NFUS = 0 if oe is not input. 


NDENS: If NDENS = 0, the densities are to vary with each 1/k. 
and are input as part of Item 14; if NDENS > 0, this 
number of densities must be input as Item 7 and each 


density value will be used for all 1/k's. 
MODES: Number of output modes, <25. 


NFREE: Number of free-free vibration modes to be input, 
< (40 - NCANT), (see next indicator). Note that NFREE 
does not vary with the number of flexible components 
(surfaces) but each free-free mode will have NSUR 
partitions and each partition will consist of the control 
point deflections for one surface. NFREE = 0 if only 


cantilever modes are used. 
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NCANT: Total number of cantilever modes to be input for 
the system. < (40 - NFREE) 


NDELM: NDELM = 0 if no rigid component generalized mass 
matrix ([Am]) is input; NDELM = 1 if [Arm] is input. 


NPUNCH: This indicator is used to obtain a printout of the dynamic 
matrix ({ U]) and to obtain the frequencies and modes in 
punched-card format. NPUNCH = 0 if no printout of 
{U] or punched output is desired; NPUNCH = 1 if only 
punched-card output is desired and NPUNCH = -1 will 
provide for the printout of [U] and the punched output. 


(The minus must be in Column 37 and the 1 in Column 40.) 


NCON: This indicator provides for changing five program test 
numbers used in the matrix iteration subroutine. 
NCON = 0 if no changes are desired and NCON # 0 if 


any of the tests are to be changed. 


3h. Data card for changes in matrix iteration tests (FORMAT 3E12.8 
and 214): There are three test numbers and two control numbers associated 
with the iteration procedure (Ref. 1, Appendix A) that are present in the 
program but that may be changed by the program user. The numbers built 
into the program are given below. An explanation of its purpose follows 


each number. 


EPSP = 0.5 1079; test for eigenvector convergence when the 
iteration procedure is approaching a single root. 

EPDP = 0.5% 107; test for convergence when the iteration 
procedure is approaching a pair of close roots. 

AITKEN = 0.9; if this test is met the program uses a procedure 


(known as Aitken's a? method) to accelerate convergence. 


oo 


NITRSP 


40; a maximum of 40 single-precision arithmetic 
iterations will be performed for each eigenvalue if its 


eigenvectors have not converged ina lesser number. 


NITRDP 


100; if the eigenvectors for any one eigenvalue have not 
converged in NITRSP single-precision iterations the 
program will then use up to a maximum of 100 double- 


precision arithmetic iterations. 


Altering any of the above numbers (NCON # 0, Item 2) requires entering all 
five numbers on one card even if some of the numbers are to be the same as 
shown above. Input EPSP, EPDP, and AIKEN in the first three E12. 8 fields 
on the card, the next two 14 fields (Columns 37 through 40 and 41 through 44) 
contain NITRSP and NITRDP, respectively. The EPSP and EPDP tests 
present in the program represent six-significant-figure accuracy although 
some loss in significance will occur as the higher modes are found; increas- 
ing these numbers would relax the requirements for accuracy and decrease 
the computing time. The AITKEN test must be less than but not too close to 
unity; an optimum convergence rate has been observed for values around 0. 9. 
The NITRSP and NITRDP can be increased to the limit of the 14 fields, and 
either control may be decreased to zero. Errors in input data often create 
convergence problems, therefore all input should be carefully checked before 
an increase is made in the maximum number of iterations. 

ITEMS 4, 5, 6, 7, 10, 13, and 14 DESCRIBED BELOW SHOULD BE OMITTED 
WHEN NAERO = 0 IS USED IN ITEM 2 (Vibration Analysis), 


4. Control card(s) for external stores and weighting matrices 
(FORMAT 1814): The consecutive fields in these cards contain ISXT), 
ISW], ISXT>, ISW2, .--> ISXTywsuR> ISWysuR: ISXT; = number of control 
points reserved for the external stores on surface 1; if no external stores 


are present on surface i, enter zero for ISXT, or leave its field blank. 
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ISW, = 0 (or blank field) if no weighting matrix is to be input for surface 
i, ISw, = 1 if a weighting matrix will be input for surface i. Include only 


the number of cards needed for i= 1, NSUR. 


5, Reference semichord, b. and 1/k (sometimes denoted by 
V/b_) series (FORMAT 6El2. 8): These reference parameters are used 
in computing the flutter velocities [Eq. (38)]. The b, (feet) may be any 
value, but the proper relationship must exist between the 1/k's used for 
reference reduced velocities and the /k 5 used when computing the AIC 
matrices for each surface; that is, V/k = (1 /k,)(b./b 3) where bo is the 
reference semichord for surface i. To input these parameters enter b. in 
Field 1 of the first card; Field 2 and consecutive fields in the first and 
successive cards (if NAERO > 5) contain the 1/k 's, < 20. 


6. The La and Ss. series (FORMAT 6E12. 8): A reference semi- 
chord (b 5) and a reference semispan (s,) must be input for the rigid 
component if NFUS = 1 and for each flexible component if NAERO > 0. 
When NFUS = 1 the consecutive fields in the first and successive (if 
NSUR > 2) card(s) contain in order Bi s_,b ex -b 


gt et Fa? rNSUR 
When NFUS = 0, the consecutive fields in the necessary card(s) 


» and 


°“NSUR’ 
contain in order boy S)> bis S5> C4 bUNSUR’ and SNSUR’ 

7. Density series (FORMAT 6E12.8): Omit this input if 
NDEN = 0 inItem 2. If NDEN > 0, begin in Field 1 of the first card and 


enter NDEN densities, <20. 


8. Rigid component generalized mass matrix [Am] (FORMAT 
6E£12.8): (Omit this matrix when NDELM = 0 in Item 2.) The size of the 
[Am] matrix for the rigid component(s) of the system must be compatible 
with the corresponding matrix for the flexible components given by 
[n, 171 M] (had; i.e. , each élement in [Am] is based upon the same rigid- 
body generalized coordinates as the respective element in the product matrix. 
If NDELM = 1, input the [Am] matrix (size NRIGID X NRIGID) by column with 


each column beginning on a new card. 
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9. Mass matrix, [MJ]: Although the [M] matrix defined in 
SectionI, Part C, has NSUR X NSUR partitions, only the nonzero partitions 
are input; i.e., a separate mass matrix ([ M],) is input for each surface. 
The sequence for considering the surfaces was established in Items 4 and 6. 


Repeat the following input for i = 1, NSUR. 
a. Size control card j (FORMAT 1814): 


Field 1 contains NSIZE., the order of [M]. (external 


stores control points + surface control points, < 49). 


b. Control card(s) for omitting zeros in the coupled mass 
matrix [M], (FORMAT 1814): 


Field 1 and consecutive fields in the first and successive 
cards (NSIZE > 9) contain in order LOW), LHIGH), 
LOW>, LHIGH?2, ---, LOWng1z7E° LHIGHystzE° meee 
is the row number in which the first nonzero appears in 
Column j, and LHIGH ; is the row number in which the last 
nonzero element appears in Column j. If only one nonzero 


element is in Column j, the row number in which it appears 


must be used for both BON: and LHIGH,. 


con Elements in [M], matrix (FORMAT 6E12.8): The elements 
are entered by column as shown below; the input for each column starts ona 


new card. 


Column 1: elements LOW, through LHIGH). 


Column 2: elements LOW, through LHIGH),. 


Column NSIZE: elements LOW, o)7p through LHIGH\c17¢- 


Any zero elements appearing in rows between sti and ery must be 
entered or their respective fields left blank. If all elements in a column 


are zero, one zero element must be entered for this column (a blank card 
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may be used) giving LOW; = sacina for this element in the above control 
cards. If external stores are present (ISXT, > 0, Item 4) all store control 
points must be considered before the surface control points; i.e., the 
elements representing the external stores mass must occupy the upper 


left-hand corner of the mass matrix. 


NOTE: One possible method for determining the [M] ma rices 
is suggested in Appendix B of Ref. l. 


10. Rigid component modal matrix, [h (Omit this input when 


Rol! 
NFUS = 0, Item 2.) The rigid-body modes for the rigid component are 

input in order to compute the incremental generalized force matrix [AQ] 
that results from aerodynamic loads on the component. The number of rows 
in [hp] must be the same as the number of control points considered when 
computing the [C.,] matrices, the number of columns must agree with 


NRIGID, Item 2. These data are input in the following manner. 
a. Size card (FORMAT 1814): 
Field 1 contains NROWS, number of rows in [h,,! : 


b. Elements in [hp] (FORMAT 6E12.8): Tabulate the matrix 


elements by column with each column beginning on a new card. 

ll. Rigid-body modal matrix, [hn]: (Omit this input if 
NRIGID = 0, Item 2.) If NRIGID > O, input the [hp] matrix by partitions; 
iné., input [h,], size (NSIZE, X NRIGID) for each surface. The [hp]; 


matrix requires the following input, repeating the input for i= 1, NSUR. 
a. Size control card i (FORMAT 1814): 


Field 1 contains NSIZE., the number of rows (<49) in 
[hp], (must agree with order of [M,] ). 
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b. [h,], matrix elements (FORMAT 6E12. 8): Enter the 


matrix elements by column with each column beginning ona new card. 


12. Vibration modes ([hy] matrix) and frequencies: The [h,.] 
matrix is also partitioned and a separate [ho]; matrix is considered for 
each surface. The columns in [ho], are the free vibration modes (orthogonal 
free-free and/or cantilever) for surface i. The total number of modes 
(NFREE + NCANT) used for the system cannct exceed 40. The frequency 


(cps) associated with each mode must also be input. 


Input the frequencies and modal column matrices ([hy] ,) in the 
following manner: if NFREE = 0 repeat Parts a, b, and c (below) for 
i=1, NSUR; if NFREE > 0 omit Part b fori = 2, NSUR. 


a. Control card A (FORMAT 1814): 


Field 1 contains NSIZE., the number of rows in [hy];, 
£49 (must agree with order of [ M],). 


Field 2 contains NMODES, number of columns in [hy 
(NFREE + NCANT,). NCANT; is the number of cantilever 


modes to be input for surface i. 


Field 3 contains NFORM; NFORM = 1 if the frequencies 
and the [hI], matrix elements are to be input using 
FORTRAN (FORMAT 6E12. 8); NFORM = 0 if these data are 


input using column binary format. 


b: Frequencies (op, series i) (Format in above control card i) 
If only cantilever modes are used for the system, an w,, series is input 
following each control card i The o's are listed in the order in which 
their respective modal columns will appear in [hy],. (Number of Wp S = 
NMODES, above card.) 


If NF REE > 0, i=1 for Op series,, that is, only one Op 
series will be input for the system and this series must follow control 


card 1 (control card for the first surface). The number (< 40) of wp's in 


= 32 


a 
ae 


the series must equal NF REE + NCANT used in Item 2. The W's are 


listed in the order in which their respective modal columns appear in 


the unpartitioned [h,] matrix. 


The Or series are input using the format specified in 


control card 2 NFORM = 1 (FORMAT 6E12. 8) permits entering six w 4's per 


card; use only the number of cards necessary to complete the beat 
NFORM = 0 (column binary format) should be used only if the W'S and 
modes are available as punched-card output from a suitable vibration 
program. When this format is used, the first Wr starts in Origin 1; the 


card(s) £2 are followed by a 


number of wp 's per card 4 22 and the Wr 


TRA card. 
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S [h.], matrix elements (Format given in control card ,): 
When only cantilever modes are used for the system, the modal column 
order used in any one [h). is independent of the order used for any other 
free-free mode is partitioned for NSUR surfaces and the appropriate part 


but if any free-free modes are used (NFREE > 0, Item 2), each 


of the mode must have the same column location in each [h.], matrix. The 
free-free modes must appear as the first NFREE columns in each [hy], 
matrix, the remaining columns (if any) are the cantilever modes peculiar 
to surface i. Use the format given in control card ; and input [hy], by 
columns. If NFORM = 1, begin each column ona new card and use the 
number of cards (49) necessary to complete the column. If NFORM = 0, 


each column begins in Origin 1 of a new card, the number of elements per 


card £22, and the card(s) £3 for each column are followed by a TRA card. 


NOTE: Punched cards output from the IBM program based on Ref. 1 may 


be used for the above Farts a, b, and c with no alterations if the output data 


—h 


“The TRA card has a 7 and a 9 punch in Column 1, Columns 2 through 72 
are blank, and Columns 73 through 80 will contain the characters used for 
identification and sequencing in the punched-card output deck. 
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(cantilever or free-free modes and frequencies) are suitable for any 


component or the complete system considered in the modal analysis. 


13. Weighting matrix, [W],, i=1, NSUR: (No provisions have 
been made for entering a [W] matrix for the rigid component; any adjust- 
ments to Pe,.o] must be made before it is input.) The presence or absence 
of a [W] matrix for surface i is determined by ISW,;, Item 4, if ISW, =1, a 
[W] matrix must be input for surface i. When the matrix is input, it 
consists of two parts if control points for external stores have been reserved 
(ISXT, > 0, Item 4). The input for [w]. matrix is as follows, with this 


exception: omit Parts a and b when ISXT, = QO. 


a Control card , for external stores [w], partition 


(FORMAT 1814): 


Field 1 contains NXST,; NXST, = 0 if no [Ww], matrix is 
input for the external stores area (the program will use a 
unit matrix, f 1); NXST, = number of control points 


reserved (= ISXT,) if the partition is input. 
Field 2 is a dummy field and may be left blank. 


Field 3 contains NFORM; NFORM = | if the [Ww], 
matrix elements will be input using FORTRAN (FORMAT 
6E12. 8); NFORM = 0 if the matrix elements will be input 


using column binary format. 


Field 4 contains NROW; NROW = 0 if the [ w] i matrix 
elements will be input by column; NROW = 1 if the matrix 


elements will be input by row. 


b. Elements in external stores [w], partition (Format given 
in control card i): Omit this input if NXST, = 0 in control card ro If 
NXST; > 0, input the matrix elements using the format combination given 


in control card 2 that is, if NFORM = 1 and NROW = 0, tabulate the elements 
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by column and begin each column ona new card, and if NFORM = 1 and 
NROW = 1 tabulate the elements by row with each row starting on a new 
card. When the column binary format is used (NFORM = 0 and NROW = 0) 
Column 1 starts in card Origin 1, Column 2 in Location (1 + IXST.), 
Column 3 in Location (1 + 2X ISXT.), etc., anda TRA card must follow the 
input for this [w], partition. 


fos Control card ; for component [Ww], (FORMAT 1814): 


Field 1 contains NSIZE., number of control points used 
on surface i. (Do not include control points for external 


stores. ) 


Field 2 contains NPART, the number of partitions in the 
component [w].. NPART = 1 for a full (unpartitioned) 
matrix. (The [w], matrix is often sparse, sometimes 
diagonal, and may be of rather large order, < 49; for this 
reason, provision is made for partitioning the matrix and 


entering only the nonzero partitions. ) 


Field 3 contains NFORM; NFORM is defined above in 


Part a. 


Field 4 contains NROW; NROW is also defined above in 
Part a. 


ns 
The input below is repeated for NPART partitions (j= 1, NPART). 
d. Partition j size (FORMAT 1814): 


Field 1 contains the order of partition j (if NPART = 1, the 


order = NSIZE, used in above control card i) 


e. Elements in partition j (Format given in the above control 


card i)! The elements in each partition must be input using the format 


“Part dis input only if NFORM = 1 (FORTRAN FORMAT) in the above 
control card; if NFORM = 0 (column binary format) the unpartitioned 
matrix is punched by columns and input as Parte. 
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combination given in Fields 3 and 4 of control card .. See the explanation 


given in Part b (above) pertaining to the use of the format combinations. 


14. Aerodynamic data: The aerodynamic input consists of NAERO 
(Item 2) sets of AICs. One set of AICs consists of the AlCs for each surface 
which have the same reference 1/k (see Item 5). If NDENS = 0 (Item 2), a 
density series will precede each set of AICs. Input the density series (only 
if Item 7 was omitted) and the AIC matrices for each surface with the follow- 


ing input order. Repeat the order for each (1/k);, j=1, NAERO. 
a. Density series for (1 /ik) 3 (Omit this input if NDENS > 0.) 
(a. 1) Control card ; (FORMAT 1814): 


Field 1 contains NRHO, number of densities to be 
entered for (1 /k);, £20. 


(a.2) Densities (FORMAT 6El12. 8): Use the number of 


cards (44) necessary to complete the series. 


b. Rigid component [c.]; matrix: (Omit this input if 
NFUS = 0, Item 2.) 


(b.1) Reference 1/k for [Cr od; (FORMAT 6E12. 8): 


Field 1 in this card normally contains the 1/k used 
in computing [G55 The 1/k card is used only to 
help identify [Co]; in the input deck and in the 
printed output; it may be a blank card, but the 
AICs in [0]; must be computed for a 1/k which 
properly relates them to the j'th 1/k.. (See Items 
5 and 6.) 


(b.2) Control card j for [o.1; (FORMAT 1814): 


Field 1 contains NSIZE, number of control points 
on the rigid component (must agree with NSIZE 
used for (hy). 
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Field 2 contains NPART, the number of partitions 


in [ey J i NPART = 1 is used for an unpartitioned 


matrix. 
Field 3 contains NFORM as defined in Item 13.a. 


Field 4 contains NROW also defined in Item 13.a. 
ote 


The input listed below must be repeated for NPART partitions. 
(b. 3) Partition k size, k=1, NPART (FORMAT 18I4): 


Field 1 in this card contains N, the order of 
partition k. If NPART = 1, N = NSIZE used in 


control card ; 


(b. 4) Elements in partition k (Format given in control 
card ha All elements in the AIC matrices are 
complex numbers; but each partition is input as 
though it is a real matrix, size NX2N. The 
imaginary parts of the elements form the 
even-numbered columns. Input each partition in 
[C1 . with the format combination specified in 
Fields 3 and 4 of the above control card. The 


format combinations are explained in Item 13.b. 


c. [c.].; matrix for flexible surface i: The input order for 
[C,] ij will differ from the order given for [o); only if external stores 
control points are reserved on surface i (ISXT, > 0, Item 4). The input 
given below must be repeated for i = 1; NSUR and Parts c.1 andc.2 are 


omitted when no external stores are present on surface i. 


= 

Part b.3 is input only if NFORM = 1 (FORTRAN FORMAT) in the above 
control card; if NFORM = 0 (column binary format) the unpartitioned AIC 
matrix is punched by columns and input as Part b.4. 


aay 


(c. 1) 


(c. 2) 


(c. 3) 
(c. 4) 


The input below must be 


(c. 5) 


(c. 6) 


“See footnote page 37. 


ee ne ~ ssa ee 


Control card ii for external stores [c,] ij 
partition (FORMAT 1814): 


Field 1 contains NXST, = 0 if no [Cc,] ij partition 


will be input. (The program sets the partition to 


zero. ) NXST, = number of control points 
reserved for external stores if the external 


stores [C,1;; is to be input. 


Field 2 is a dummy field and may be left blank. 
Field 3 contains NFORM as defined in Item 13.a. 


Field 4 contains NROW as defined in Item 13.a. 


Elements in external stores [Cc.],; partition 
(Format in above control card i; The elements 
in the partition are complex numbers and will be 
input in the same manner as the elements in the 
Leia! partitions (Part b.4, above), with the 
exception of format which is given in the above 


control card. 


Reference 1/k for surface i (FORMAT 6E12. 8): 
(See Part b.1, above.) 


Control card ,, for surface [on]; partition 


(FORMAT 1814): (See Part b.2, above.) 
repeated for NPART partitions.” 


Partition k size, k = 1, NPART (FORMAT 1814): 
(See Part b.3, above). 


Elements in partition k: (Format in above control 


card, Part c. 4). 
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Note: Punched cards output from IBM programs based on the sources listed 
below may be used (with no alterations) for Item 14, Parts b.1, b.2, b.3, 
and b.4 or Parts c.3, c.4, c.5, and c.6 when the theory is appropriate to 
the problem. Punched cards output from the program based on slender -body 
theory may be used for Item 14, Parts c.1 and c.2 if the first card (1/k.) 

is removed from the deck. The list shown is complete only at this writing; 


additional compatible sources are currently under development. 


W. P. Rodden, E. F. Farkas, H. Maicom, and A. M. Kliszewski. 
“Aerodynamic Influence Coefficients from Incompressible Strip Theory: 
Analytical Development and Computational Procedure."" Aerospace 
Corporation Report No. TDR-169(3230-11)TN-5, 3 September 1962. 


W. P. Rodden, E. F. Farkas, H. Malcom, and A. M. Kliszewski. 
‘Aerodynamic Influence Coefficients from Supersonic Strip Theory: 
Analytical Development and Computational Procedure.'' Aerospace 
Corporation Report No. TDR-169(3230-11)TN-1, 1 August 1962. 


W. P. Rodden, E. F. Farkas, H. Malcom, and A. M. Kliszewski. 
‘Aerodynamic Influence Coefficients from Piston Theory: Analytical 
Development and Computational Procedure.'' Aerospace Corporation 
Report No. TDR-169(3230-11)TN-2, 15 August 1962. 


W. P. Rodden, E. F. Farkas, and G. Y. Takata. ''Aerodynamic 
Influence Coefficients from Slender-Body Theory: Analytical Development 
and Computational Procedure. '' Aerospace Corporation Report No. 


TDR-169(3230-11)TN-6, 31 October 1962. 
Cc. Example Keypunch Forms 


Example keypunch forms are given on the following pages. Columns 
73 through 80 are reserved for data deck identification and sequencing. Only 
the cards with sequencing are used in the two data decks set up for the 
example problems; the lines with Columns 73 through 80 left blank are for 


description of input data. 


39 = 
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SECTION IV 


PROGRAM OUTPUT 


A. Printed Output 
be All input data. 
25 The generalized mass matrices (iid or [mie]) corresponding 


to the input free vibration modal matrices (igs or [ht] yc 


3. The dynamic matrix if NPUNCH = -1 is used in the general control 
card. 

4. For the vibration analysis or for each 1/k in the flutter analysis: 

a. The eigenvalue for each output mode followed by the number 


of single- and double-precision iterations and the number of Aitken 


accelerations. 


b. The eigenvectors for each mode followed by the check 


eigenvalues and eigenvectors. 


c. The frequencies in cycles per second and, ina flutter 
analysis, the structural damping coefficient and the velocity (knots) associ- 


ated with each frequency. 


d. The final flutter or vibration modes for each component 


normalized on the largest element in each mode. 


e. If NPUNCH = +1, the sequencing numbers used in 
identifying the punched cards containing the frequencies and the final modes 


will appear below the respective printed data. 


Di In a vibration analysis, the generalized mass corresponding to 
each output (free vibration) mode will follow the above printed output [see 
Eqs. (40) and (62)]. 


6. A number of different statements may be printed indicating 
machine- or program-detected errors in input data such as wrong format or 


incompatibility in the number of rows in [hp] ? [M] ? and [nh] i 
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7. If the program or the machine fails in the iteration portion of the 
program, a note will be printed stating the type or cause of failure. Some of 
the various program stops and machine failures are explained in Ref. 1, 


Appendix A. 


8. If convergence is not obtained in the allowable number of itera- 
tions, a note will be printed and the program will continue. (In this case, 
the eigenvalues and eigenvectors should be compared with the check eigen- 
values and eigenvectors to determine if the convergence is sufficiently 


accurate.) 


9. The printed output for the two example problems is shown on the 


following pages. 
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00 JOEETILED*°O- 
00-349 T86E 12 °0- 
00-3STIZISLE°O- 
00 3894968¢6°0 

10 366666660°0O 

00-386219802°0 

00-J6E6EOEZE “0 

00 391969LL9°0 

00-3129084%0T°0 

6 NWN103 

10 366666660°0 

00 3€66€6986°O 

00 35298%1%9°0 

00 36494198€9°0 

00-3TO2TILL2°0 

00-JEZ20S8EL2°0 

20-J66E6898E°O 

Z20-F9ILEBEBIE °0 

00-3%S62szTt°0- 
00-384S0860T °0- 


€ NWN109 


S3dOW Y311NI4 NIG 


20-398€0LL68 °0- 
°0 
00-319€068TT°O- 
10-302599669°0- 
00-3SL040€E9T°0- 
00-3SL244LET°0- 
10-32419998S °0- 
10-36€922618°0- 
ZO-3SETISYEE*O- 
10~F09SLE9BZ °0- 
8 NwN109 
10-380%9026T°0 
°o- 
10-3$1924299 °0 
10-3L992991S °0O 
10-3S22S08¢6°0 
10-318269948°0 
10-3990416S6°0 
10-392982Z106°0 
10-30621TE296°0 
10-398T10826°0 


2 NWN109 


JIVAYNS 


00 316680€%6°0 
10 366666660°0 
TO-JLTL192L28 °0- 
00-3LETELOZT°O- 
00 ATSLZTOTS°0- 
00 3IBOLEEIES*O- 
10-3Z19L922L°0- 
10-31S5066219°0- 
TO-J6ETZTELE°O 
T0-3124S2SSS°0 
L NWN1109 
oo 3TIS9OTLE°O 
10 366666660°0 
00 360STLS4S°0 
00 32S92SELS°0 
0O-J*STESBET*O 
OO-JETLOBE9T*O 
00-J1SZOvTET °O- 
00-395S9680T °0- 
00-38292ZE92°0- 
00-301221922°0- 


a NWN109 


ANANTH OM- DAO 


ANMSENOM OAS 


- 


onl 


af 


40 32v6T09ZE°O G06-368694892°0- 10 JLLE0SS56°0 S 

40 J0LSZELEZ°O 00 302800219 °0- 10 3981S58669°0 % 

90 368SZ29TST°O 00-326T6T8L2°0 10 3260€TL44°0 € 

€0 J012@L4626°0 10-32S22902T °0- 10 3568260%12°0 2Z 

€0 36906 T002°0 00 JITI9E9%E9°0- TO 3S8629902°0 T 
({S1ONH) ALIDOWA ONT dwvcd (S43) v9SWO 300W 
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B. Punched Output 


le A deck of punched cards is output from the vibration analysis 
portion of this program that is suitable as input to other programs requiring 
free-free vibration modes and frequencies. The flutter modes and frequencies 
also may be punched in cards if it is so desired. If a printout of the dynamic 
matrix is needed (NPUNCH = -1) but no punched-card output is wanted, the 


machine operator must be instructed to ignore the output punch tape. 


2. All punched output is sequenced in order on Column 73 through 
80 starting with HM140000. Each mode will have NSUR (number of surfaces) 


partitions, and the output data will be in the following order. 
a. Control card for surface 1 (FORMAT 1814): 


Field 1 contains NSIZE number of control points on 


1’? 
component l. 

Field 2 contains NMODES, number of output modes 
punched. 


b. Frequencies (column binary format): The frequencies are 
punched in the order in which they appear in the frequency column in the 
printed output. The column binary card(s) (only one card if NMODES < 22) 


containing the frequencies are followed by a TRA card. 


cor Modes for component 1 (column binary format): The 
component 1 partition of each vibration mode is punched as a column with a 
TRA card following the output for each column. If the final flutter modes 
are punched in cards, a TRA card will follow both the output for the real 


column and the output for the imaginary column forming each mode. 


d. Control cards and modes for component i, i= 2, NSUR: 


Parts a and c (above) are repeated for NSUR surfaces when NSUR > 1. 
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SECTION V 


PROCESSING INFORMATION 


A. Operation 
FORT RAN II MONITOR system (uses chain job capability). 
B. Estimated Machine Time 


Because of the many variables involved such as number of input modes, 
number of surfaces, number of iterations required for eigenvector conver- 
gence, etc., it is not possible to give a practical metnod of estimating 
machine time. Complete but independent decks may be stacked to save 
program read-in time. The two data decks for the example problems were 


stacked and the machine time was 1.84 minutes. 
Cc. Machine Components Used 


The following machine logical tape units were used; each is set equal 


to the symbolic tape unique to the program. 


TAPE 2 (input tape) = NTAPE2 
TAPE 3 (output print tape) = NTAPE3 
TAPE 11 (chain tape) = NTAPE1 
TAPE 12 (output punch tape) = NTAPEO 


TAPES 13, 5, 6, 7, 8, and 14 (utility tapes) = NTAPES4, 5, 6, 7, 8, 


and 9, respectively. 


NFEAPE1 contains two chains. Chain 1 uses six utility tapes for storing 
matrices; the number of rows and number of columns precedes each matrix. 
The symbolic utility tapes and the matrices stored on each are listed below 


(i = 1, NSUR on all tapes except NTAPE 9). 
NTAPE4: 1) [M], 


If NRIGID >0 
2) [h,] for the system 


=7B< 


If NRIGID > 0 and NFREE > 0 


3) fa,» 
NTAPES: 1) Pagel: 
NTAPE6: If NRIGID > 0 and NFUS > 0 
1) [hpol 
If any ISW; >0 
2) [w], 
NTAPE?7: 1) [w].[c,]. or [c,]; 
NTAPESB: If NRIGID > 0 


1) [M], = (tu, + 32. 1749 ,b¢.s,[w]{C,]; ) 
j = 1, NDENS or 1, NRHO 
2) (MI ,lhp], 
3) {h] [M4], 
NTAPE9: 1) (myyl = ln, [MI [hp], . i= 1 for [hy], 
if any free-free modes are input; 
If NRIGID > 0 
2) F_e), = [bg]; [Ml,{hy, 
Mp; = [nply [M] [hp]; 
If NRIGID > 0 and NFREE > 0 
4) [amllap el 
5) lappl [aml 
6) [appl lamllagsyl 
If NRIGID > 0 
7) [Fy_l 7) = ((agl TMI [ng] + fam] )"? 
4): [its eg)? 
Chain2 also uses NTAPES8 for storing computations usedinthe MITERS subrou- 


tine and for storing the normalized vibration or flutter modes for each surface. 
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SECTION VI 


PROGRAM NOTES 


i. Subroutine Used 


In addition to the main program (page 91), the following subroutines 


were used. 
1. RDLN, reads and prints the title card (page 127). 


2. MMULTD, computes the product [C] of two matrices [A] 
and [B] (page 129). 


3. MPRINT, prints matrix in matrix format (page 132). 


4. MREAD, reads a matrix in either column binary or FORTRAN 
format (page 135). 


5. BINRD, reads column binary cards from tape (page 140). 

6. MNVRSX, computes the inverse of a complex matrix (page 150). 
7. INVERS, computes the inverse of a real matrix (page 153). 

8. BINPU, writes column binary cards on tape (page 158). 

ele MPUNCH, punches column binary cards (page 168). 

10. MITERS, real or complex matrix iteration (page 172). 

seid NPNRMX, vector normalization (page 185). 

12. SWEEPX, computes the true mode and sweeps it from the 


related matrix (page 190). 


13. CLOSE, matrix iteration using double-precision arithmetic 
(page 196). 
14. DPFORM, converts double-precision numbers to single- 


precision numbers and conversely (page 207). 


15. DPMLTX, matrix multiplication using single- or double - 


precision numbers (page 211). 


LT 


B. Generalized Tapes 

The Aerospace IBM 7090 computer logical tape units used in the 
program are TAPES 2, 3, 5, 6, 7, 8, 11, 12, 13 and 14, but these may be 
altered by placing the desired machine logical tape units on symbolic cards 
HM140716 through HM140725. These symbolic cards appear in the main 
program, and the logical tape units are defined in the comments preceding 


the cards. 


ae. 


FLOW DIAG 


=,¢ © 
OMENSIONS 


F Py NAERO = IF \oo 
EQUIVALENCE HRs NAERO-! NAERO 
- FORMATS 


NSRFLX = 
NFREE + 
NCANT 

NSURFS= 

NSUR + NFUS 


lm: ReSxST,, SW ROOREF, | 
a \=|,NSUR Fe, MAE 


NPOINT= 
ZNSIZE ISUR 


SECTION VII 


FLOW DIAGRAM -- MAIN PROGRAM 


(og 


DEFINE PROGRAM 
"AND CALL 


NAERO = DEFINE RDUN 
NAERO-| TAPES ITERATION 
NTAPE6 CONSTANTS (TITLE CARD) 


NSRFLX = : . 
PRINT: NFREE + gel ae Re 
1ON NCANT z= i 
CORRTANTS NSURFS= CANTILEVER FREE-FREE 
NSUR + NFUS MODES MODES 


TERATION 
ONSTANTS 


A READ: NC=2 READ: 


© 


F ah 
NAERO 


READ: 
ISXST, ,ISW BREF, I/Kr; R2=2* NRIGID BRj Sj 
iy i i= 1, NAERO i NR2 RIG] 12 |, SURFS 


1=1, NSUR 
° 


PRINT: 
TITLE- 
MASS 
MATRIX 


; REWIND 
‘ * He es DEFINE, PROGRAM CALL READ: 
_ DEF! NTAPES icone ROLN CONTROL CARD 

_ TAPES NTAPEG Gunconties (TITLE CARD) 


PRINT: 
Thee = 
FREE-FREE 

MODES 


NC=2 
NR2=2*NRIGID 


Si 
NSURFS 


tl, 


PRINT: 
TITLE- 
MASS 

MATRIX 


KI=KL 
+NSIZE 


PRINT: 
4 vibe) READ: 
TITLE- NMO CNFREE,NSRFLX, 
rew=(0] FLEXIBLE eo 0 
MODES 


NMODES = 
NSRFLX 
K2=1 


CALL MMULTD 


c= (ne) [M] [he] 


HT = HT + ARF 
FFM=FFM+C 


CALL MMULTD 


ARF = [*nturmce, 
Prey, EMI) 


nsarcx] A 


© 


{ sprinc} = 
[NAGONAL Mer] 


2 NO i] hl tn 
wuss | -@ es) 


—\+ O 


CALL MMULTD 
-l 
G=-[Am] *ARF 


CALL INVERS 
va: (Am) 


CHAIN | (CONT) 


READ: 
[e;. weade | 


READ: 
wi. 
iz |, NSRFLX 


CALL MMULTD 
save= [M][hp] 


{spring} 2 
[DIAGONAL Mer] 


HT = HT + ARF 
FFM=FFM+C 


CALL MMULTO 
c= [ne }"[M] [he] 


- {sprinc} = 
“DIAGONAL Mrr] 


CALL MMULTD 
= 
G= -[am] alii 
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Se, 


a 


PRINT: 


00 


NO 
178 
WEIGHTING 
1=1,NSUR MATRIX 


REWIND 


STORE: 
NTAPE6 ('78) NTAPES6 
NTAPET 


Cw] | NTAPES 


READ: 7 READ: 
NXST, IFORM, EXTERNAL | 
IROW STORES 


dQ=[0] 


(READ: 
Poa) 


READ: 
NXST,IFORM, 
IROW 


WEIGHTING 


STORES 
MATRIX 


© ~ @ @ 


REWIND 

NTAPEG6 
NTAPE? 
NTAPES 


NRHO = NVEL = 
NDENS NVEL +1 


D Ce READ oO CALL MMULTD | 
| 
[soees"y [cn], F=[ChILhng] 


STORE: READ: STORE: 
NTAPE7 NTAPE6 NTAPE7 
tw], Mw, [c], 


CHAIN | (CONT.) 


READ: -READ: REWIND: 
NXST,IFORM EXTERNAL AP 
IROW , [Stones | i ns 


a 9 
RHO , NVEL = -.+ 
+ |,NRHO NVEL +! i= 


CALL MMULTD 


F=[Ch] ER, ] 


READ: CALL MMULTD STORE: 
NTAPES hate NTAPE7 
(w; i Mi, [c], 


- © ~~) » & 
{ 
bo CON= RHO (IRHO) + DRRAR = CONDNDG i 
300 174 
TRHO=/, NRHO HB eS HSE BARMRR = DMBAR 


i 
) =) 9 (z=) (xs) | 
Ce Ose 1 ¢ 


STORE: 
NTAPES 


Ure)" [ln] 


STORE: REAO: CALL STORE: 
Pracl LA] [oy [a] r=[re}"[8) [hau] [>e} [i [omg] 


BARMRR rea = 
BARMRR 
BARMRR wae2S 


K=ISUR+NFUS 
CON=RHO (IRHO) 
3 ES, 32.174 


REWIND: 
NTAPE4, 
NTAPES, 
NTAPE 7, 
NTAPES 


CON=RHD (IRHO) 


DMBAR = CON #00 
3B, * HS, 332.176 OM 


+ 
BARMRR -OMBAR 


save =[0] 


STORE 
NTAPES 


[m)" [a] 


CALL 
MMULTD 


A= [ma] LM] [om] 


CALL TORE 
MMULTD NTAPES 


SAVE=[onr|” a=[one]" HT [onr]” [Om] [one] 


YES 


Vn Est cen mn OTTO Am OO 


STORE 
NTAPES 


[ Msoun] [Maa] 


S) 2 BARMRR aca = oo 
BARMRR 
BARMRR wae: 


266 
TSUR=},NSUR 


YES 
aa a 
YES 
ves 


CHAIN (1, 11) 


K«TSUR+NFUS 
CON =RHO (IRHO} 
HF 965, M3274 


i (w) _ & 
«[rv]"[m) 


o 


“i wll aon game. La "ie 


[one]" [Ai] [+n] 


ERRO! 
INVERSE 


ee J > vES 

BO ALL STORE: iets ‘iis 

266 MMULTD x 2) NTAPES___ (268) (00) ai 

Loi) *LR od) \\ [nn Me : 
YES 
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wee 


MODES > MODE 
DIMENSIONS REWIND wae: [0] 
EQUIVALENCES NTAPES #FF:[o} 


COMMON FORMATS ner: fo} 


HFA, HER, ARF 


WRF, =< HRE, +aaF 
Le L +NCensze 


Ao KONE oO NSIZE 


STOME [Merl = PRINT 
rare? (i/[sPrins TITLE 
(rye PREGITIF-M ourpur oars 


cat 
OPMLT x 


LH «HEA « HAR 


CALL 
pPMuT x 


HRRA* HRA & HAE 


[Mr] 


DAMP? Ac / hy 


M4 


TI,NCe MOOES, VEL* PICON way « 


bret # ¥/ba 


TILL NENA 
NORMAL IZE [ae]8¥ LaAGesT 
ELM (US) IN EACH MOD 


INDEX=(SUR MiQAD= H¢L.L2 US 
usrisur), L+us HIWOEX LI M¢2,L1« 


2 USUNDEX +14.) 
(SUR +t LT 


CALL DPMLTK CALL DPMLTX 


' 
rr eal 
#2 Be HF 2 


HF slit eplo[n] 


MODES = MODE 


> REWIND: Do REWIND : HRF = [0] REAO: 
ore tS NTAPE? ° 450 NTAPES neF=[o] NTARES 
MAXR= 49 IRHO = 1, NRHO, arfe-|[mM 
NueTe) HER = [0.] FFI sun 


READ: 


READ: 
HFF = HEF +ARF NTAPES HRF=HRF +ARF 


ARF = [Pee], 


HRF, = HRF, +ARF 
L=L+NC-#NSIZE 


HFRy= HFRy + ARF 
K=K+NC *NSIZE 


NTAPES 
ARF= [one], 


sree meena emererenteeererrstir 


iF PRINT: S1ORE 
PRINT STORE: 
GALL 0 (1/[SPRING TITLE - S OYNAMIC NAKSR=0 CALL MITERS wrore’? 
Ff DPMLT x (1) FREQ (TI) oureur pate MATRIX [y] NAKDR=O 1, MODES te] 
os 2 i! 
[Fae] ANF =HFR » HRRR ( [mer] -[arr]) 


PRINT 
TITLE «f iWaL 
PLUTTER 
“ODES 


READ: 
NTAPES 


HF! +[rr] a 


DAMP © Ai Shy 


VEL=PICON @ wy * 
Peat * ¥/be 


COLL OP MTX 
Hes Hey «[h] 


378 
ISU «1 ,NSUR 


ves Mii, US LIU) 


NTAPES HILLISUSTIA} WOKDR IL) ISUR 


nO, 
CALL NPNRMX NORMALIZE 
fr JOY MeL US TIBOR) 


ISUR= |, MoU 


al 
+ 
CALL OPMLTK pursa-ane ye? uS+ U5 +HF 
He Behe = 


a CHAIN 2 


HFR=HFR+ ARF 


MA WINd, 
iarena 


READ, 
WIAPEN 


rae 


mrt 


CALA GiratT® 
nr em al) 


mie HE 


we [>.] 


ie 


A 
0, Pe VEC 


< BUH ed. ) 
teu ‘ 
vsti use.) dan eumieune!..3"}) 
A 


oo melibie 
see 


(SERi 2 NIA, 


|ves 


|. NPAMMER NORMAL PE (3) 
jBY #O.dyustisuaa) : 


d __| 


PRINT? 
NGO=98765 
. GENERALIZED 
GALL DPMLTX MAXR = 50 


USZHFF er MASSES, — US; 


a 


Ss 


SECTION VIII 


SYMBOLIC LISTING 


Some of the program FORTRAN symbols which were not defined in 
Section III, Part B are defined below. 


FORTRAN Symbols 
BARMRRi(I, J) 
BOOLT 

BREF 

BR(I) 

DM(I, J) 
DMBAR(I, J) 
DQ(I, J) 


MAXS, MAXR, MAXP, 
MAXQ 


NFLEX, NFLEX2 
NGO 


NPOINT 


NRIGID, NR2 
NSIZES(I) 
NSRFLEX 
NSURFS 
PICON 


RHO(I) 


Definitions 
—_— — -1 
[mpp! and isc 
BCD titles 
b, (reference) 


b.. series 
ri 


[4m] 


[Am] 
gal [ Chol [hp] 


Dimensional number of rows in working 


matrices 
Total number of flexible modes 
Chain 1 entry control number 


Total number of control points on all 


flexible surfaces 

Total number of rigid-body modes 
Number of control points on surface i 
NFREE + NCANT 

NSUR + NFUS, (NFUS = 0 or 1) 

2x 7x 0.5921 


p., i= 1,NDENS 
1 
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SPRING 


T, G, A, SAVE, H, ARF, 
C, FFM, F, HR, HT 


VELCT Y(I) 


The complete symbolic listing 


[K/4 x? | 


Working areas of core 


(1/k.), = (V/bw). 


is given on the following pages. 


90 - 


T d9Vd 


SS90%1WH 
+S90%1WH 
ES904THH 
2S904%1KH 
TS90%THH 
0S90%1WH 
64900 1WH 
849051WH 
L79041WH 
9*90%1WH 
S*90%1THH 
+e90%1WH 
EY90"THH 
Ze90"THH 
T9904 THH 
0%9041WH 
6E90%1WH 
SE90%1nH 
LESO4TWH 
SESO%1 HH 
SE9O04TWH 
E904 1WH 
EESO4HTHH 
ZESOSTHH 
TEVO+IWH 
O£90%TWH 
6Z90%1WH 
8Z9041WH 
L290471WH 
9290%1WH 
6290%1hKH 
9Z90%1WH 
EZ9041HH 
ZZ90%1WH 
TZ9041WH 
OZ9041THH 
6190%THH 


€9/80/0T 


*(AOULING (LZISILITD S(CdSULING (9TSILT) '( dda (STSI1T} 
"(dSd3'(VISILID S(SGUVINS (ETSIET) S(O3dVIN'O(ZTSILT) 
#CTHANS (TIS) LIDS (ALITIAS (G69) 11) SCONTUAS £(0S9)1T) 
®(NODId* (644}1T) S(OHUNS(B4F)LT) S(GuaVAS(Lo4) 17) 
*(HUWUVE (SLEILT) SCuvend *(eOEd1T) ‘(bd 
S(TEZILI) S(WOS(SOTILID SCSSCHLTILT)D SC uet(ESTIIT) 
*(XTSASNE(Z2STILID §(S3NASCTSTILID S(B3usS(TITILT) 
*(QHYS (TOIL) 4(1100G*(6B8)1T) SCHINNdN' (88011) FINI IWAINGS 
(OONS (28211) S CONS (98P1T) SCZ2XFTAN' OSS )ET) 
SCZUNS (FSET) S(SSUMSNECEBILID SCdxVW (28 d1T) 
S(SXWWECTBDLI) SCOxvWS(OBDLT) SCuxvWS(62)11) 
*CXDTINS (SL)LT) 8 CANIGANS(LZZ)LT) S(S300W'( 92017) 
*(SNSON'(G2211) S(SNANS (HL)LT) SCANVON‘S (ELIIT) 
®(QAMANS(ZLPLT) SCGTOTUNS(TLILI)D SCunsN{(OL)1T1) 
®(63dVIN'(O9)1T) S(83dVINGE(89PLT) S(L3dVINEO (LOLI) 
S(9FdVIN'(99DLT) S(G3dVENS(S9DLT) SC 43dVIN'S(99)1T) 
*CEIdVINGCEIILT) S(Z3dVINGCZ9DLT) SC TIGVIN'CTSIITDS 
(SIZISNS(THDLT) SCOTIMSTECTZILTD SCCTIESKXST S(T ILI) «INS IWAINOS 


ana tr NOPE 


aANmMT NOP OD 


(0SOS)W4d* (OZIALITSAS (0540S) 3AVS '(O4IONI US * 

(OOTSOS)L S$ (0449) 4uv S(ZTSOSIH S(O0TS9)1H 

#1054999 *(O0T'OS)S* (94GG0TIYH *(GS‘0S)9 ‘(GGT ‘ds)¥ 

*(OZG)4T  §{02)S3ZISN *(OSIHSIHT *(0S)M01 

*(OHIOSYs SC Z2TE9PUMWUVE S(ZTI9)YVEWNG §(2T'9)0G 
*(o*9)wa *(OZIOHY *(TZIS 8(1Z)Us FCAZIMST *(0Z)1SXST NOISNAWIG 


aN 


COTOLTaKad 
*GIOIUNIOO “(G1OLYUN#Z*GIOTUN) YVGWG ‘(GIOTYN‘GISTYNIWG 
* (SJGOWN#UNSNIOINA S(CISTUNSSZISN#UNSN) HH 
*(3ZISNIHOIHT S(3ZISNIMOTSCOISTUNIWO SC T+YNSKIS 
*©(T4UNSNIYS “(SNIGNIGHY “CUNSNIMST *CUNSNILISXST NGISN3IWIG 
*39V4YUNS JNO ANV NO LAdNI S3GOh YISWNN IW101=SITOHWN 
G3MO011V 39V4xNS JNO ANV NO SINIOd TOWLNCD Y3AWNAN WI01=3Z1SN 
GaMoTW S31G08 GIO9IY 30 YIAWNK IWLOL=GI9T YN 
*G3MONIV SSILISN3G JO Y3EWNN TVLOL=SNION 
*G3M011V S39VSUNS 430 YIGWON IWLOL=HNSN 


ANOS 


"SIV ONISN *GOHLAW IWGOW V AG SISATWNVY YSLINTS ONY NOLIVUGTA 


YVYUOUVUYUYYUYUOY 


=9i1.= 


€6904TWH ONIAHOISM ON HET *X9 SHIT *Z2IT 33VdUNSHB *XBE HT } i1VWYOs TZ 
26904%TWH ( OBZLVNIWYSL WVYOONd *3NILNOW 3SYSANI NI YOUN THL%) 1VWYOS 02 


T690%TWH ( XIYLVW SSVW LNSNOdWCD OISTY H6Z *XO€ OHT) LvWYOd 61 
0690%1TWH ( *SINIOd WOWLNOD HLT 1 
6890%1TWH *‘OllT *S300W LNSNOdWOD OIDIY HEZ *XOE OHT) AVWHOs BI 
88904 1WH (€ 8°023T = ¥4 w/°T HOT T 
1890%1WH *xs *x€ XIMLVW DIWVNAGOUZY HOZ ‘XEE THT} =ivwuod LT 
9890%1TWH ( SLNIOd HZ I 
G89041WH TONINOD H8*6IT *XIVLVW O¥3V LNINOdWOD OIDIY H6Z ‘X82 THT) AVWYCd OT 
4890%TWH { XIYLVW ONTLHOISM H8T *X9% THT) 1VWYO3 ST 
€890%TWH (| XIMLVW TWOOW 3T7@IX3173 HEZ2 *X€% THT) Lvwx03 41 
Z89041WH (SINIDd TOMLNOD HST “911 “HT $211 39V4dUNSHB ‘x8E€ HT) LvWYyOd €l 
18904%TWH ( XITuLVW TVOOW AGOG AISIY H¥Z*xz% THT) LvWuod ZT 
0890%TWH ( "°**O3NNILNOD WVYDD%d $35N9V LON OO (HLE 4 
6LOOFIWHS4IT } $09199dX3 SEINIOd TOYLNAD JO Y3SEWAN TVLOL ONV (H8Y I 
82904TWH SoIl )} *XIMLVW SIHL SLNIOd JO¥LNOD 3O A3AWAN OHT4) 1VWHOs TT 
2190%1TWH ( XITULVW SSVW H2T *x8% THT) ivWudd OT 
GLOOFIWH( (8°02 31 *XS% HI) // S3LTINFJINDIYS HET “9IT *XS% OHT) IwWdO4d 6 
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